Garbay- Toulouse Phase Transition in Heisenberg Spin-Glasses 

in Three Dimensions 
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We examine three-dimensional ±J Heisenberg models with and without random anisotropies 
in a magnetic field. We calculate both the stiffness exponent 6 S at absolute zero temperature 
and spin-glass correlation lengths for the longitudinal and transverse spin components at finite 
temperatures. We suggest that, contrary to a chirality scenario predicted by Kawamura and his 
co-workers, a Gabay- Toulouse phase transition occurs when the anisotropies are absent, although 
no phase transition occurs when they are present. 
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Phase transitions of three-dimensional (3D) Heisen- 
berg spin-glass (SG) models have attracted much interest 
in recent years. Two phase-transition scenarios are in- 
volved in controversy. One is the SG scenario, in which 
a usual SG phase transition takes place at a finite tem- 
perature. However, this scenario has been believed to 
hold only when anisotropies are present [l[ . The other is 
a chirality- glass (CG) scenario, proposed by Kawamura 
and his co-workers 0, H, H|- In the CG scenario, not the 
spins, but the local chiralitics freeze at a finite temper- 
ature. In this scenario, the SG phase transition never 
occurs in isotropic SG models. The freezing of the spins 
was suggested to occur through coupling of the spins and 
the local chiralities by random anisotropies. 

Kawamura and his co-workers gave three pieces of ev- 
idence of the CG scenario in the isotropic case: the 
stiffness exponent for the chiralities is positive X > 0, 
whereas that for the spins is negative S < 0;2]; only 
the chirality autocorrelation exhibits a pronounced ag- 
ing effect at low temperatures Q; the chirality overlap 
distribution P(q x ) exhibits a one-step-like replica sym- 
metry breaking (RSB) behaviorQ. However, reexamina- 
tions of those properties revealed different aspects: the 
stiffness exponent for the spins in a lattice with open 
boundaries is positive S > || H, Q; the spin auto- 
correlation of a system, in which a uniform rotation is 
removed, exhibits an agingeffect similar to that of the 
chirality autocorrelation § and the SG susceptibility 
Xsg exhibits a divergence behavior to a finite temper- 
ature TsgQ- Recently, the SG phase transition tem- 
perature Tsg and the CG phase transition temperature 
Tcg were estimated using a non-equilibrium relaxation 
method [ldl and the scaling method of the SG correlation 
length [lH. Results suggested that Tsg = 2cg in both 
methods. Based on those results, the SG scenario has 
come back also in the isotropic models. That is, a usual 
SG phase transition takes place at a finite temperature. 
Freezing of the local chiralities results. 

The controversy surrounding those two scenarios has 
reached a new stage. It has been speculated that true 
SG properties are visible only in large lattices, e.g., the 
L x L x L lattice with L > 20, because the coupling of 



the spins and the local chiralities, which exists even for 
the isotropic model in small lattices, loosens for L —> 
oo[l2j. Campos et al. quite recently studied the model 
for big lattices (L — 24 and 32) [l^ to resolve this issue. 
Their results suggest that the lower critical dimension 
di of this model is equal to or a slightly smaller than 3 
(dl ;$ 3) and that a large finite size correction exists in 
the scaling property of the model with d = 3. Having 
taken into account this correction, they also suggested 
that Tsg = Tcg- However, objections exist in relation to 
their interpretation[3|. Unfortunately, it is too difficult 
to resolve this issue herein. 

Two scenarios predict different aspects for a finite mag- 
netic field H ^ 0. In the SG scenario, a usual phase 
transition will take place, which is characterized by a 
freezing of the transverse component of the spin, i.e., a 
Garbay- Toulouse (GT) phase transition[l^]. In the CG 
scenario, the CG phase transition will occur, but the SG 
phase transition is absent 17f . More interesting is a case 
in which anisotropies are present. In the SG scenario, 
the SG phase transition will disappear because of a ran- 
dom field effect [l6;]. On the other hand, Imagawa and 
Kawamura predicted that the CG transition still occurs 
at H ^ 0, accompanied with the SG phase transition[l8j. 

In this letter, we present an examination of the phase 
transition of the ±J Heisenberg models with and with- 
out random anisotropies at a finite magnetic field H =/= 0. 
Special attention is devoted to an induced magnetic mo- 
ment (Si) at each site i. We consider a SG spin compo- 
nent, Si(= Si — (Si)), to examine cooperative phenom- 
ena of the system. Results show that, in the isotropic 
model, the ground state stiffness and the scaling prop- 
erty of the SG correlation length suggest the presence of 
the SG (GT) phase transition, like those at H = 0. On 
the other hand, in the anisotropic model, both the CG 
transition and the SG transition disappear. Therefore, 
we suggest that a usual SG phase transition occurs in 
the isotropic Heisenberg SG model at H ^ as well as 
at H = 0. 

We study the ±J Heisenberg SG models in three di- 
mensions (d = 3) in a magnetic field H described using 
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the Hamiltonian: 

w = - E JijS i s j -EE Drs?s? -hJ2s: 

(ij) (ij) liu i 



(1) 



where Si is the classical vector spin of \Si\ = 1; Jij = + J 



or — J with the same probability of 1 /2. The second term 
expresses the anisotropic energy; D^"{= Dj" = D^) 
{(J,,f = x,y,z) are symmetric random anisotropic con- 
stants distributed in the range [-D : D], The lattice is a 
simple cubic lattice of L x L x (L+l)(= AT) with periodic 
or skew boundary conditions along two L directions and 
a periodic boundary condition along the (L + l) direction. 
We consider two cases: (A) D = and (B) D + 1 0. 

(A) Isotropic case of D = 

First, we consider the ground state stiffness of the 
model using a method proposed by Matsubara et al.\jj^. 
Here we consider lattices of L x L x (L + l) with open 
boundaries for the (L + 1) direction. The lattice has two 
opposite surfaces Oi and fix,+i- We first determine the 
ground state spin configuration {Si = sj + S^-} and its 
energy E\. Then, fixing all the spins on the surface fii, 
all the spins on the surface are rotated by the same 

angle <f> = tt/2 around the z-axis and fixed. Under this 
boundary condition, we calculate the minimum energy 
of the system, E^, which is always higher than E\. The 
stiffness of the system might be characterized by the ex- 
cess energy AEl(= E^ — E®). The stiffness exponent 9 S 
might be defined by the relation AE L oc L 6s . 

We have calculated [AE^] of the model up to 
L = 12, together with the parallel (the magnetiza- 
tion) and the transverse components of the spins, [M z (= 
| £\ sf |/AT)] and \S A -(= £\ \S^\/N)], having used a ge- 
netic algorithm 19], in which [• • • ] means a sample aver- 
age. Numbers of samples with different bond distribu- 
tions are 100 ~ 1000. Figure 1 shows [M z ] and [S^] as 
functions of H. Those values depend little on L, sug- 
gesting that they are those for L — > oo. In fact, [M z ] 
exhibits a characteristic property of the SG, i.e., it in- 
creases rapidly with H and saturates gradually at high 
magnetic fields H s ~ 7 J. Consequently, [S- 1 ] has a con- 
siderable value up to H s . Figure 2 shows [AEt] for 
several H in a Log-Log form. Using least-squares fit- 
ting, we estimated the stiffness exponent as 8 S = 0.58 ± 
0.02, 0.65±0.02, 0.63±0.03, 0.63±0.06, and -0.18±0.15, 
respectively, for H/J = 0,2, 4, 6, and 7. It is interest- 
ing that 6 S are positive and almost equal for H < H s . 
This result is analogous to that in the spin-flop (SP) 
phase of an antiferromagnetic Heisenberg (AFH) model, 
in which S = 1 for H < H s (= H c = 12 J). Therefore, 
we expect that the SG (GT) phase transition occurs at 
H < H c (— 6.5 ± 0.5J). 

We next examine the phase transition of the model. 
We consider a two-replica system with {Sf} and {Sf }. 
At H ^ 0, the spins are polarized to the ^-direction: 
TOj(= (S z ) a — (S*)P) + 1 0, where (• ■ • ) means a thermal 
average. Magnitudes of will vary from site to site. 




FIG. 1: Magnetization [M z ] and the transverse component 
[S ± ] for various lattice with L as functions of H. 
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FIG. 2: (Color online) Excess energy [A El] for various lat- 
tices with L. The symbols denote, from the above, those at 
H/J = 0,2,4,6 and 7. 



Figure 3 shows their distribution P(m,i). In fact, dis- 
tribute in a very wide range. The slight size dependence 
of P(rrii) reveals that rrij is a purely magnetic-field in- 
duced one. Then we subtract m 4 (= (0,0, m*)) from the 
original spin [111]: Sf' 13 = Sf'^ 3 — rrii. Hereafter, we call 
S"' 13 SG components and consider their cooperative phe- 
nomena. 

We consider the spin-glass correlation lengths for 
the longitudinal (rj = ||) and transverse (77 =-L) compo- 
nents. We calculate them using a formula |llj: 
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2sin(fc min /2) \XsG( k 



x n SG (Q) 



- 1 



1/2 



(2) 



where k min = (0, 0, 2iz/{L + 1)). The fc-dependent 
SG susceptibility is given as x| G (k) = N [(\<j zz (k)\) 2 ] , 
and x^ G (fe) = Nj2^ =x y [(|^(fc)|) 2 ], with q^(k) = 

W J2i S^S^ exp (ik ■ Ri). If a SG phase transition oc- 
curs, the correlation length divided by the system size L, 
has the following scaling property: 



fV . 



T SG (H)) 



(3) 



where v is the correlation length exponent, Tsq(H) is the 
transition temperature at H, and ^ represents a scaling 
function. 
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FIG. 3: (Color online) Distribution of the site magnetizations 
riii of the model with D — in a finite magnetic field. T = 
0.25 J is slightly higher than T SG (~ 0.2 J) at H = 0. The 
arrow indicates the average value of mi for L = 13. 



We performed a simulation of this two replica system 
on the lattice with L < 23 having used a temperature 
exchange MC method [2(| with an over-relaxation [l3j . 
Numbers of samples are 128 for the largest lattice; equi- 
libration is checked by monitoring the stability of the 
results against runs at least twice as long. Figures 4(a) 
and 4(b) respectively show ij^/L and ^\/L at H = 0.2 J 
as functions of T. These two quantities exhibit differ- 
ent size dependence. £^ / L for different L cross around 
T = 0.22J, suggesting the presence of the phase tran- 
sition. In fact, choosing T§q ~ 0.215J, we can scale 
£l (see Fig. 5) using a finite size correction exponent 
of tj> = 0.9 predicted by Campos et a/.[l3[. This result 
is compatible with the ground state study. In contrast, 
(} L / L seem not to cross, even at low temperatures, sug- 
gesting the absence of the phase transition in the longi- 
tudinal component. This result is also compatible with 
that of the Ising SG model at H ^ 0, where no-AT line 
was suggested [21j. 

We have considered properties of the isotropic SG 
model at if ^ 0. Both the ground state stiffness and the 
scaling properties of the SG correlation length suggest 
the occurrence of the SG phase transition in the trans- 
verse component of the spins. That is, the GT phase 
transition will occur at H ^ 0. 

(B) Anisotropic case of D ^ 

Next we consider the anisotropic model with D = 0.1J. 
We first note that the system has only the spin reversal 
symmetry of {Si} — > {—/£,} at H = 0. This symmetry 
breaks at H ^ 0. Consequently, the system at H ^ 
will exhibit no phase transition associated with the spins. 
On the other hand, Imagawa and Kawamura (IK)|18| 
examined the CG phase transition and suggested that the 
CG phase transition occurs because of the one-step-like 
RSB accompanied with the SG phase transition. They 
considered cooperative phenomena of the original spins 
{Sf } and {Sf }. We have reexamined it^22| using SG 
components {Sf} and {Sf}. Figure 6 shows ^/L as 
functions of T. In stark contrast to the case of D = 0, 
the £l/L for different L seem not intersect at any finite 
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FIG. 4: (Color online) Plots of (a) £f /L and (b) Q/L of the 
model with D — 0. 
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FIG. 5: (Color online) A finite-size scaling plot of fj". 



temperature. This result supports the former argument. 

Does the CG phase transition really occur at H ^ 0? 
We have also reexamined the chirality transition using 
the SG components {Sf} and {Sf}. Fi gur e 7 shows 
the chirality overlap distribution of P((j , x )[23(] at a low 
temperature. In marked contrast to the IK results (see 
Fig. 8 in [18]), P{q x ) exhibits a single peak at q x — 
0, which becomes sharper as L increases. This result 
indicates no freezing of the local chiralities. 

In summary, we have examined the phase transition 
of the three-dimensional ± J Heisenberg models at finite 
magnetic fields H ^ 0. When anisotropies are absent, re- 
sults suggest the occurrence of the SG (GT) phase tran- 
sition as in those at H = 0. On the other hand, no 
phase transition occurs when they are present. These re- 
sults support the SG scenario of the phase transition of 
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FIG. 6: (Color online) Plots of £jf / L of an anisotropic ± J 
model with D = 0.1J at H = 0.2J. 



Heisenberg SG models. 

Finally, we should note that present results seem also 
to be compatible with experimental observations. Pe- 
tit et al. performed torque experiments of a series of 
Heisenberg SGs with different magnitudes of local ran- 
dom anisotropics in magnetic fields -ff[24|]. They found 
an irreversibility of the transverse spin components below 
a finite temperature Ti(H), i.e., a putative SG transition 
temperature. They showed that, when the magnitude of 



the anisotropy is very weak, Tj (H) is almost independent 
of H at low magnetic fields. It is suppressed strongly 
at H ^ as the magnetude is increased. If Ti{H) is 
a crossover temperature between the paramagnetic-like 
state and the SG-like state[25j, our present results could 
explain these observations. 




FIG. 7: (Color online) The chiral-overlap distribution func- 
tions at T = 0.16 J for different lattice sizes L. 
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